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ABSTRACT 

In this Letter, we consider the role played by the sedimentation of helium nuclei on the emis- 
sivity and metallicity distribution of the X-ray emitting plasma in the cores of relaxed galaxy 
clusters. We model the gas density and temperature profiles of nearby cooling core clusters to 
estimate the gravitational acceleration acting on the helium and show that its sedimentation 
time scale is too long with respect to the present age of these objects to play a significant role. 
However, we argue that these time scales have to be definitely lower in the past allowing the 
helium to settle in the cluster cooling cores and raise its abundance to values higher than the 
solar one. A direct consequence of this speculation is that the helium, by increasing the total 
X-ray emissivity, reduces the measured metal abundance in the inner (r < 20 kpc) cluster 
regions. 
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1 INTRODUCTION 

The plasma in X-ray emitting clusters of galaxies is composed by 
hydrogen (H), helium (He) and various ionized metals of which 
iron (Fe) is the most interesting from the astrophysical point of 
view owing to its connection to the star formation activity (e.g. 
Renzini 2003). In an X-ray emitting plasma with solar abundance, 
for each atom of hydrogen a number density of 9.77 x 10^^ ions 
of He and 4.68 x 10^^ ions of Fe is expected (as in Anders & 
Grevesse 1989; for the most recent determinations in Grevesse & 
Sauval 1998 and Asplund, Grevesse & Sauval 2005, 8.51 x 10"^ 
ions of He correspond to each atom of H, whereas the number den- 
sity of Fe is 3.16 x 10"^ and 2.82 x 10"^, respectively). Being 
the atomic weight A of (H, He, Fe) = (1,4, 56), the mass fraction 
on the total is 0.75 for the hydrogen (generally labeled as X), 0.24 
for helium (called Y) and the rest (Z = 1 - X - y 0.017) for 
heavier elements. 

Diffusion of helium and other metals can occur in the core 
regions of the intracluster plasma under the attractive action of 
the central gravitational potential, enhancing their abundances on 
time scales comparable to the cluster age. Qin & Wu (2000) and 
Chuzhoy & Nusser (2003), following the work by Fabian & Pringle 
(1977), Rephaeh (1978), Abramopoulos et al. (1981), Gilfanov & 
Syunyaev (1984) on the sedimentation of elements in clusters, stud- 
ied the concentration of helium in cores concluding that it can grav- 
itationally settle within a Hubble time. Abramopoulos et al. (1981) 
obtained the same results of Qin & Wu assuming a global thermal 
equilibrium and the Boltzmann equation, while Gilfanov & Syun- 
yaev (1984) applied the diffusion equation to the light elements 



(atomic number Zi < 16) where the changes are more efficient. 
Therefore, in a H-He plasma, the helium tends to be more centrally 
peaked than the hydrogen. 

More recently, Chuzhoy & Loeb (2004) solved the full dif- 
fusion equations derived by Burgers (1969) for a multicomponent 
fluid, finding that the diffusion, if the suppression due to magnetic 
fields is modest, can alter X-ray properties as the steepening inward 
of the baryon distribution as well as the spectrum and evolution of 
stars forming out from helium-rich gas. However, all these studies 
assume a gas temperature constant in time and space in evaluat- 
ing the gravitational acceleration in the cluster core, where instead 
steep spatial gradients are observed in cool core clusters. 

In the present work, we analyze the effects of sedimented el- 
ements. He in particular, on the description of the X-ray emitting 
plasma that is generally assumed with a fixed solar abundance of 
helium (as tabulated in, e.g., Anders & Grevesse 1989; on more 
recent estimates of the solar chemical composition see, for exam- 
ple, Asplund, Grevesse & Sauval 2005) and a relative contribution 
from metals determined through the equivalent width of the de- 
tected emission lines in X-ray spectra. We model properly the cen- 
tral steepening of the observed gas temperature and density pro- 
files to estimate the gravitation potential in the inner cluster regions 
where the effect of the gravitational sedimentation is expected to be 
more relevant. In these regions, we prove that the underestimation 
of the amount of helium can propagate to inaccurate measurements 
of the total X-ray emissivity, gas density and metallicity. 
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Figure 2. (Left panel) Equipartition time between H and He (thick dotted line) and between H and Fe (thin dotted line) in the Centaurus cluster; sedimentation 
times, tged = / dr /vgcii{r), of He (solid line). The dots mark the radius rj^ at which the material is accumulated from the regions beyond with a 
sedimentation time represented by the line. (Right panel) Sedimentation times of He for the three objects under examination (from thickest to thinnest line: 
Centaurus, A2199, A1795). The cluster potentials are described by eauationl?! given the best-fit values in Table I. A flat temperature profile is assumed (see 
text at the end of Section 2). The shaded region ranges between 2 X lO'' years, which is about the cooling time associated with a gas with solar abundance of 
helium and proton density of 0.2 cm~'^ at ~ 1 keV. and the age of the Universe for the assumed cosmology, 1.3 X 10^'' years. 



2 SEDIMENTATION OF HELIUM 

Spitzer (1956) presents the physics of fully ionized gases that we 
adopt in the following discussion. In general, the time of equipar- 
tition between H and He (Fe) ions is low when it is compared 
to the cluster age, whereas the sedimentation time is compara- 
ble to it. Hereafter, we assume a cosmology (Ifo, f^m, J^a) = 
(70 km s~^Mpc~^, 0.3, 0.7) that implies an age of the Universe 
of 1.3 X 10^° years. 

To describe the action of the gravitational acceleration in drift- 
ing inward the plasma elements, we consider an intracluster gas 
in hydrostatic equilibrium with a Navarro, Frenk & White (1995, 
hereafter NEW) potential 



with X = r/rs. We are interested to measure such potential in the 
inner regions of a cluster, where a steep gradient in temperature is 
generally observed. A model for the gas density profile obtained 
analytically from the hydrostatic equilibrium equation (HEE) with 
a NEW potential is available in literature (e.g. Makino, Sasaki & 
Suto 1998, Ettori & Eabian 1999), but it assumes an isothermal gas. 
We relax this assumption to accommodate the observed gradients 
and modify accordingly the model of the gas density profile. We 
assume a temperature profile with the functional form 
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that is just the rearranged expression of the equation that Allen, 
Schmidt & Fabian (2001) show to be able to reproduce the tem- 
perature profiles of relaxed clusters (see Vikhlinin et al. 2005). We 
use then this model in the HEE with a NEW potential to recover a 
numerical expression of the gas density profile: 



t{x) 



dx 



(3) 



ngas(K) = ngas(O)-^ exp -r? 

with T] = AnGpsrsfimp/Tn^^^. 

These functional forms of the temperature and density profiles 
are then joint-fitted with a minimization to the deprojected data 
of Centaurus (Sanders et al. 2002), A2199 (Johnstone et al. 2002) 
and A1795 (Ettori et al. 2002). The best-fit results are presented in 
Table 1 and plotted in FigQ 

A proper fit of the gas density and temperature gradient allows 
us to define accurately the cluster gravitational acceleration g by 
using the NEW potential in equation^ 



-fir) 
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with the parameters rj, Tmax and Vs — r/x obtained as best-fit re- 
sults of the deprojected gas density and temperature profiles. By 
using, instead, the best-fit parameters from an isothermal NEW gas 
density profile, as generally done for this kind of analysis, one 
would infer a potential in the inner regions that is larger than the 
adopted value by a factor of ~ 4, mainly owing to the assumed 
higher temperature value that propagates linearly to the gravita- 
tional acceleration. 

Eor a Boltzmann distribution of particles labeled 1 with den- 
sity ni and thermal velocity ?;th ~ (2fcT/^imp)^^^ in a plasma 
with temperature kT in hydrostatic equilibrium with a NEW poten- 
tial g{r), the drift velocity of the heavier ions 2 with respect to 1 is 
given by (e.g. Spitzer 1956) 

3ml AiA2 < g{r) 



cd(r) = 
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where log A is the Coulomb logarithm equals to 37.9 + 
In ((fer/10keV)(ne,i/10"^cm-^)-i/2). Typical velocities of 
about 2 kpc Gyr~^ are expected which can make reasonable the 
sedimentation of helium in the inner 10 kpc or so. The time elapsed 
in drifting from rinn and rout is the sedimentation time, t^cd = 
jj-out t;r/iiscd(r), that is therefore proportional to ni kT~^^^. 

The time required to sediment helium in the inner cluster 
regions is in the order of few Giga-years (see also Chuzhoy & 
Loeb 2004). Owing to the competing drag forces from the pro- 
tons moving upwards and the helium diffusing inwards, heavier 
elements can sediment with similar velocity, as discussed in Chu- 
zoy & Nasser (2003). We conclude that sedimentation of helium is 
not an efficient process, given the present condition of the cluster 
plasma. 

However, the intracluster medium evolved to a well-defined 
cool core during its formation history. We can assume that the gas 
started originally with a flat temperature distribution equal to the 
value measured in the outskirts and a gas density profile obtained 
from the equation|3|by fixing t{x) = 1 and by requiring that the 
outer pressure remains unchanged as the overall gravitational po- 
tential. While the evolution of a cool core from an initially isother- 
mal plasma is the natural consequence of the radiative processes 
taking place in the central regions at higher density, there is not 



Table 1. Best-fit results of the gas density and temperature profile models 
as plotted in Fig. HI 





Centaurus 


A2199 


A 1795 


z 


0.0104 


0.0309 


0.0632 


Ta (kpc) 


16 


33 


246 


ne(0) (cm-3) 


0.264 


0.133 


0.089 


To (keV) 


0.48 


1.36 


2.25 


Tmax (keV) 


3.92 


5.22 


12.24 


b 


0.91 


0.94 


0.86 


V 


12.1 


7.6 


15.9 


X^(ne) 


212 


12 


29 


xHt) 


453 


20 


10 


X^/dof 


665/10 


32/10 


39/12 



yet detailed studies from, e.g., cosmological simulations of struc- 
ture formation that discuss the evolution of cool cores owing to the 
still unknown balance of energetic losses and feedbacks that rules 
the X-ray emitting plasma in the cluster cores (for example, the 
temperature profiles are expected to be flatter and density higher 
at 2 = 1 with respect to the local measurements in Ettori et al. 
2004, but Bums et al. 2004 suggest an alternative scenario in which 
cool cores in massive clusters are formed through the merging of 
subclumps with its own cool core). These assumptions imply that 
the gas temperature increases in the inner regions by a factor of 
~ 3 and the density decreases in the center by about an order of 
magnitude. As shown in Eig. |5| these changes in gas temperature 
and density values cause the sedimentation time to be shorter. This 
makes reasonable the scenario in which (i) the helium sedimented 
during the initial collapse and subsequent relaxation of the clus- 
ter plasma and (ii) it resides nowadays in super-solar abundance in 
cluster cores. 



3 EFFECTS OF THE SEDIMENTED HELIUM 

Intracluster plasma emits mainly by X-rays through bremsstrahlung 
with a typical emissivity. 



e (X A(Tg, 



)nc ^ MtZfn,, 



(6) 



where A(Tgas) is the cooling function of a plasma with tempera- 
ture Tgas, Zi is the electric charge (f.e. 1 for H, 2 for He, 24 for 
Ee), Mi is the element abundance relative to the solar value and 
Hi is the density of the species in exam. It is thus evident that 
any metal enhancement of its relative abundance with respect to 
the hydrogen, as occurs when sedimentation is relevant, affect the 
total X-ray emission. It is important to note that the X-ray emis- 
sivity is generally measured under the assumption that the helium 
abundance is fixed to the solar value relatively to hydrogen, i.e. 
e oc 1.85np. In the extreme case of a plasma composed by pure 
helium (X — 0,Y = 1), no = 2nHo and e oc Stv^^, with the 
direct implication that if this plasma is assumed with solar abun- 
dance, the total density is overestimated by a factor of 2. Moreover, 
the change in the relative abundance of H and He induces variation 
in the number of electrons and ions present in the plasma 
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Figure 3. Deprojected metallicity profiles. An inverted gradient is present in tlie inner cluster regions. 



where Mi is the element metallicity and Ni is the particle number 
density according to Anders & Grevesse (1989), and affects the 
measurements of the mean molecular weight, 
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that is generally fixed to a value of ~ 0.6 appropriate for solar 
abundance. While cm and /i are quite independent of the metal 
component of the plasma (e.g. /i — 0.618 and cm- — 1.209 for a 
solar abundance, /i — 0.613 and cm — 1.200 for a metallicity of 
O.SZq for the reference values in Anders and Grevesse; note that, 
on the basis of the compilation in Grevesse & Sauval, jj, = 0.605 
and 0.600 and cm = 1.182 and 1.174, for a solar and 0.3 times 
solar abundance, respectively), their values are highly sensitive to 
the variation in the helium abundance: for A/hc of, e.g., 2 and 10, 
= 0.699 and 1, respectively, cm = 1.405 and 2.968. 

Through these quantities, several other cluster measured prop- 
erties are affected, like, e.g., the cooling time tcooi ~ ngTg/e ~ 
(1 + l/cM)rgas/(npAcooi), the gas mass density pg = ^(nc + 
np)mp = + CM)npmp, the total gravitational mass measured 
through the hydrostatic equilibrium equation, Aftot <x fi~^, the gas 
mass fraction, /gas = A4as /Aftot oc fj,^{l + cm)- The gas frac- 
tion, in particular, increases by a factor from 1.4 to 4.8 for 2 and 
10 times the solar abundance of helium, respectively. To quantify 
the effect of sedimented helium on, e.g., the cooling time we have 
first used the model vmekal in XSPEC (Amaud 1996) to estimate 
the cooling function Acooi ~ e/(nonp) for a range of gas tempera- 
tures, fcTgas, and He abundances, Muc, in the interval 1 and 10 with 
respect to the solar chemical compilation in Anders & Grevesse 
(1989; see command abund in XSPEC: uhc = 9.77 x lO^^nn). 
As expected, higher helium abundance induces larger emission and 
shorter cooling time of the emitting plasma, with a cooling func- 
tion that raises for, e.g., a gas at 2 keV and proton density of 0.2 
cm"^ by a factor of 1.2 (Afec = 2) and 2.6 (A/ho = 10) and 
a cooling time that decreases from 7 x 10^ years to 5.6 x 10^ 
(Mhc = 2) and 2 x 10^ (A/hc = 10) years. Under the assump- 
tion that bremsstrahlung dominates, tscd/tcooi oc ri^^^/T^^^ which 
makes evident how the formation of a cool core, with the central 
gas that becomes progressively higher in density and lower in tem- 
perature, reduces significantly any process of sedimentation. 

Finally, we have simulated typical thermal spectra of long- 
time exposed cluster cores with the ACIS-S configuration on the 



Chandra X-ray satellite generally adopted to fully exploit the larger 
collecting area at low energies with respect to ACIS-I. We have 
used the vmekal model in XSPEC to fake spectra with helium 
abundance in the range 1-10 times the solar abundance as in Anders 
& Grevesse (1989). Then, we have fitted to the simulated data both 
mekal and vmekal models with fixed solar abundance of helium, 
leaving only the normalization, the temperature and the metal (one 
single parameter in mekal, 9 parameters -O, Mg, Al, Si, S, Ar, 
Ca, Fe, Ni- in vmekal) abundance free to vary. Our results, shown 
in Fig.|4| indicate how an unrecognized super-solar abundance of 
helium causes underestimates of the metal (iron) abundance and 
overestimates of the model normalization owing to the induced in- 
crease of the total bremsstrahlung emissivity (cf. equation |6j. We 
do not measure any relevant variation in the gas temperature and 
do not find any significant dependence of these trends upon the as- 
sumed values of column density, normalization and temperature. 
An important implication of these results is that it can potentially 
explain the drops in metallicity observed through X-ray spectra ex- 
tracted from the inner (< 20 kpc) regions of several nearby bright 
galaxy clusters (e.g. Centaurus, Sanders et al. 2002; A2199, John- 
stone et al. 2002; A3581, Johnstone et al. 2005; A2204, Sanders 
et al. 2005; note that also A1795 shows a flat, marginally decreas- 
ing metal abundance in the inner two bins as shown in Ettori et al. 
2002; see few examples of this inverted gradient in Fig.|3}. Figure|4| 
shows that an enhancement of a factor of few (~ 2 — 4) in the he- 
lium abundance with respect to the assumed solar value in Anders 
& Grevesse (1989) is sufficient to explain the observed reduction 
of 20-50 per cent in the metallicity measurements. 



4 DISCUSSION AND CONCLUSIONS 

We have studied the effects of the sedimentation of helium nuclei 
on the X-ray properties of galaxy clusters. To this purpose, we have 
estimated the gravitational acceleration by adopting functional 
forms for gas density and temperature profiles that have been 
obtained under the assumption that the plasma is in the hydrostatic 
equilibrium with a NEW potential. These functional forms are 
fitted to the observed deprojected gas density and temperature 
profiles and the best-fit results are used to model properly in the 
inner cluster regions the gravitational acceleration. The observed 
gas density and temperature values do not allow the metal to settle 
down in cluster cores over timescales shorter than few 10® years. 
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Figure 4. Best-fit spectral results obtained by fitting mekal and vmekal XSPEC models to simulated Chandra spectra. The panels show how the metal 
(iron) abundance, the model normalization and temperature change as function of the helium abundance considered here between 1 (the default value) and 10 
times the solar metal composition as in Anders & Grevesse (1989). The results are independent of the input values of the absorbing column density, spectra 
normalization and metal abundance. 



On the other hand, by assuming that in the same potential the gas 
that is now describing a cool core was initially isothermal, the 
sedimentation times are reduced by 1-2 order of magnitude within 
0.5r-s ~ 0.2r2oo- The sedimentation of helium can then take 
place in cluster cores. Even modest enhancement in the helium 
abundance affects (i) the relative number of electron and ions 
CMi, (ii) the atomic mean molecular weight /i and all the X-ray 
quantities that depend on these values, such as the emissivity 
e, the gas mass density pg oc fj,{l + cm), the total gravitating 
mass, Mtot oc fi~^, the gas mass fraction, /gas oc fj?{l + cm)- 
Moreover, we show that if we model a super-solar abundance 
of helium with the solar value, we underestimate the metal 
(iron) abundance and overestimate the model normalization or 
emission measure. For example, with A/hc ~ 3, the measured iron 
abundance and emission measure are ~ 0.65 and 1.5 times the 
input values, respectively. Therefore, increased helium abundance 
might explain the drops in metallicity observed in the core of some 
nearby systems. Other implications of the suspected and plausible 
super-solar abundance of helium in the cluster cores are: (i) the 
helium-rich gas will then cool out and form helium-rich stars that 
have a typical lifetime shorter and burn hotter than Sun-like objects 
by an order of magnitude (Chris Tout, priv. comm.); (ii) the central 
elliptical galaxies of cool core clusters might be then overabundant 
in helium (see speculation on this expected characteristic in 
Lynden-Bell 1967), (iii) the central infall of helium induces 
heating of the surrounding plasma by ^ mucnucVscdg ~ 6.5 x 

IO-^^Mhc ( „ , ) ( „, g % , ) ( , erg/s/cm=^ 

VO.Olcm^^/ \ 10^"crg/g/cm J \ Ikpc/Gyr / / 

that can offset the central cooling of about nuncA{T) ~ 

2 X 10"^*^ ( oicm-a ) erg/s/cm'' by a negligible amount. On 
the other hand, the diffusion of helium can be suppressed by 
the action of confinement due to reasonable magnetic fields, as 
recently pointed out from Chuzhoy & Loeb (2004), or limited to 
the very central region (r < 20 kpc) by the turbulent motion of the 
plasma. 
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